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I. INTRODUCTION 



Effective field theory (EFT) provides a powerful framework to study low-energy phenom- 
ena in a model-independent way by exploiting the separation of scales in physical systems 
Q, [|. Only low-energy (or long-range) degrees of freedom are included explicitly, with the 
short-range physics parametrized in terms of the most general local (contact) interactions. 
Using renormalization, the influence of high-energy states on low-energy observables is cap- 
tured in a small number of constants. Thus, the EFT describes universal low-energy physics 
independent of detailed assumptions about the high-energy dynamics. Recent applications 
of EFT methods to nuclear physics have made steady progress in the two- and three-nucleon 
sectors ||, ^ [7|, In this paper, we continue a complementary program to apply EFT 
methods to many-fermion systems, with the ultimate goal of describing nuclei and nuclear 
matter. 

We showed the promise of using effective field theory for the nuclear many-body problem 
in Ref. ||, where a dilute, uniform gas of fermions interacting via short-range interactions 
(which could be highly singular potentials such as hard cores) was analyzed. In general, this 
problem is nonperturbative in the potential, so the conventional diagrammatic treatment 
sums particle-particle ladder diagrams to all orders and thereby replaces the bare interaction 
by a K matrix [[U], p| . Each K matrix is in turn replaced by an effective range expansion 
in momentum, and in the end one finds a perturbative expansion of the energy per particle 
[see Eq. (0)] in terms of the Fermi momentum kp times the S'-wave scattering length a s 
(and other effective range parameters in higher orders). 

In contrast, the EFT approach is more direct, transparent, and systematic. The EFT 
automatically recasts the problem in the form of a perturbative Fermi momentum expan- 
sion. The freedom to use different regulators and renormalization schemes was exploited 
by choosing dimensional regularization with minimal subtraction (DR/MS), for which the 
dependence on dimensional scales factors cleanly into kp dependence solely from the loop 
integrals and dependence on effective range parameters solely in the coefficients. This al- 
lowed power counting by simple dimensional analysis and manifested the universal nature 
of the low-order contributions. 

Power counting refers to a procedure that identifies the contributions to a given order in 
an EFT expansion; this is an essential feature of the EFT approach. In the dilute case, the 
power counting with DR/MS is particularly simple: each diagram contributes at a single 
order in kp and there are a finite number of diagrams at each order. This renormalization 
scheme also simplified a renormalization group analysis to identify logarithmic divergences 
that lead to nonanalytic (logarithmic) terms in the kp expansion of the energy ||, 12]. 

Thus, the application of EFT methods to the dilute Fermi system exhibits a consistent 
organization of many-body corrections, with reliable error estimates, and insight into the 
analytic structure of observables. EFT provides a model-independent description of finite- 
density observables in terms of parameters that can be fixed from scattering in the vacuum 
or from a subset of finite-density properties. The universal nature of the expansion is also 
a key feature; any underlying potential, probed at long distance, is reproduced by the same 
form and the differences lie only in the low-energy constants. This perturbative analysis 
is directly applicable to systems of trapped atoms and provides a controlled theoretical 



laboratory for studying issues such as the status of occupation numbers as observables (13 
or the nature of the Coester line HT3 



However, bound nuclear systems and the most interesting phenomena in atomic systems 
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require a nonperturbative EFT treatment. The most conspicuous features of nuclear inter- 
actions leading to nonperturbative physics are the large S"-wave scattering lengths. The 
analysis in Ref. || assumes that the (spin-averaged) scattering length is of natural size (e.g., 
of order the range of the underlying interaction). When the scattering length is large, the 
expansion applies only at very small kp, and an alternative power counting is needed to 
extend the EFT to larger densities. This power counting, developed in Refs. JTH] and ITS 



for two-particle scattering in free space, prescribes that leading order for large a s must sum 
all diagrams with non- derivative four-fermion (Co) vertices (although in higher orders the 
other vertices appear only perturbatively). For two-particle scattering, this infinite set of 
diagrams is easily summed as a geometric series. In three-body systems, large scattering 
lengths necessitate three-body input at leading order to remove regularization dependence 



1% fiq] . The shallowness of the deuteron and triton bound states ensures that the physics 
of large scattering lengths is dominant in these systems. 

At finite density, the class of diagrams with Cq vertices is much larger than in free space, 
since it also includes tadpoles, particle-hole rings, hole-hole rings, and so on. These diagrams 
are not simply summed and have not been calculated numerically. As an alternative to a 
direct numerical solution, we can seek additional expansions. Possibilities include geometric 

19| , strong-coupling, or large- N expansions. In the present work, we investigate the most 



immediate extension that is nonperturbative in kpa s by adopting an effective action formal- 
ism, which is a natural framework for implementing nonperturbative resummations. As we 
illustrate below, the loop expansion in the auxiliary field formulation for spin-independent 
forces is equivalent to a systematic 1/N expansion, where the relevant "JV" is the spin-isospin 
degeneracy g. 

In Coleman's classic lecture on the 1/N expansion, he identifies two reasons to pursue 
such expansions [p0| . First, they can be used to analyze model field theories so that intu- 
ition beyond perturbation theory can be developed. They provide nontrivial, but tractable 
examples (at least qualitatively) to build intuition about phenomena such as asymptotic 
freedom, dynamical symmetry breaking, dimensional transmutation, and nonperturbative 
confinement. We seek similar controlled insight into nonperturbative effects in nonrelativis- 
tic many-body systems. The large-g expansion has a well-defined power counting that sums 
certain classes of diagrams to all orders to provide a systematic expansion for systems where 
kpa s is small (natural scattering length) but gkpa s /7c is order unity or greater. The expan- 
sion for natural scattering length has the interesting feature that it is nonperturbative in the 
medium while the description of free-space scattering is perturbative. For large a s , the ex- 
pansion may not be fully systematic, but could still provide insight into the nonperturbative 
many-body physics. 

The second reason for pursuing 1/N expansions is that they may be relevant to physi- 
cal systems of interest. Coleman (and many others) considered quantum chromodynamics 
(QCD) with a large number N c of colors, arguing that the large N c world is in many ways 
close to the real world of N c = 3. But while 1/N expansions of QCD and other relativistic 
field theories are common |^T|, they have been exploited much less widely in nonrelativis- 



tic many-fermion physics and very little in applications to nuclear systems. In fact, the 
principle applications have been to one-dimensional Fermi systems [23|]. However, there 
have been some applications of 1/N expansions in relativistic approaches to nuclear sys- 
tems. In Ref. |[24|| , e.g., the Walecka model was studied in a 1/N expansion by extending the 
SU (2) isospin symmetry to SU (N). To our knowledge, the present work is the first on effec- 
tive action EFT and the 1/N expansion for nonrelativistic Fermi systems with short-range 
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interactions in three dimensions. 

In Ref. [E3 , expansions in 1/g are discussed, but with the caveat that they are not likely 



relevant to nuclear systems, even though g = 4 would seem to be large enough. This is 
based on the observation that with empirical nonrelativistic interactions, exchange (Fock) 
contributions are larger than direct (Hartree) contributions, whereas the former should be 
suppressed by 1/g relative to the latter. However, this dismissal may be premature. In 
a different phenomenological representation of the problem, using covariant interactions, 
the Hartree pieces do, in fact, dominate ||2b]| . Furthermore, the expansion may be useful 



to describe part of the physics, such as if the long-distance pion physics is removed E7 
or averaged out f28fl . Indeed, analyses of phenomenological energy functionals fit to bulk 
nuclear properties suggest a robust power counting with short-distance scale A roughly 
600 MeV [p9|j . At nuclear matter equilibrium densities, this would imply kp/A is less than 
one, but gkp/A is greater than one! (However, the expansion parameter may actually be 
g hp /An, as suggested by the analysis in Sect. |TV| .) 

An interesting special case of the g — > oo limit follows if we also take the Fermi momentum 
hp to zero, with the density (which is proportional to gk F ) held constant. We refer to this 
procedure as the "Bose limit," because it generates the ground state of a dilute Bose system 
(under the assumption that the ground state evolves adiabatically from the noninteracting 
state). This limit was noted long ago in Refs. |3(J] and |31| but was not exploited until 



Brandow used it in comparing fermionic ( He) and bosonic ( He) many-body systems |32] 



In Ref. p5 |, Jackson and Wettig used the Bose limit in analyzing minimal many-body 
approximations and cleanly derived the leading corrections to the Bose energy (see also 
Ref. in this connection.) 

In Sect. II, we summarize the effective action formalism for fermions using an auxiliary 
field. In Sect. PTT| , we review the EFT treatment of a dilute Fermi system with short-range, 
spin-independent interactions and then apply the effective action formalism. We carry out 
the loop expansion, which is seen to correspond to a 1/g expansion. In Sect. [TV], we calculate 
the energy per particle explicitly to next-to-leading order (NLO) and perform a stability 
analysis of the ground state. The Bose limit is considered in Sect. [V|, and Sect. |V| contains 
a summary of our results and conclusions, along with plans for further investigations. 



II. EFFECTIVE ACTION FOR FERMIONS WITH AUXILIARY FIELDS 



In this section, we review the effective action formalism for fermions with auxiliary fields 
at zero temperature (T = 0). We adopt the notation and spirit of the general discussion 
of effective actions in Ref. EM and largely follow the specific treatment of nonrelativistic 
fermions by Fukuda et al. [33]. Our discussion will be somewhat schematic in order to 
focus on the new aspects of the EFT treatment. Details of the renormalization of the 
effective action and caveats related to convergence of the path integrals, Wick rotations 
from Euclidean space, zero-temperature limits, and so on are well documented and can be 
found in the standard literature |20|, ^4[ [36|, [37j . 

Consider a system of fermions in an external potential v (x) interacting via a spin- 
independent, local two-body interaction Uq(x — y). (The extension to many-body forces 
in the context of an EFT is outlined in Appendix [A].) Such a system is described by the 
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Lagrangian 



C tpl(x) I + ^ + ^ + v(x) ) ip a (x) - ^ipl(x)il;^(y)Uo(x - y)ip p (y)i) a (x) 



where a, j3 are spin indices, fi is the chemical potential, and M the fermion mass. We define 
a generating functional Z[J] and an energy functional E[J] via the relation (we follow the 
notation from Ref. W- 



Z[J) = e - iE[J] = J ThfiTtye 



(2) 



where J(x) is an external source. 1 For simplicity, normalization factors are considered to be 
implicit in the functional integration measure. 

The strategy is first to reduce Z[J] to a Gaussian integral in the fermion Grassmann 
fields by "integrating in" an auxiliary field Using the identity [R5] 



JVaexp (§ fd 4 x Jd 4 y \a(x) - ipl(x)ip a (x))U {x - y)[a{y) - ^l(y)ipp(y)] 
J Vaexp(^Jd 4 x fd*ya(x)U (x-y)cr(y)) 



(3) 



we introduce an auxiliary field o with bosonic quantum numbers and obtain for the gener- 
ating functional with J = (i.e., the partition function) 



Z[J}\ 



V ip^Vip Va exp i I d^xip^x) (id t + + fi + v{x) ^ip a {x) 
-il>U*)Mx) fd 4 zU (x-z)a(z) + ^a(x) I ' d 4 y U (x - y)a{y)] , (4) 



where the denominator of Eq. has been absorbed into the normalization of Z. The 
Vifj^T>tp path integral is now Gaussian and leads to a determinant of the operator between 
if)' and tp, which we identify as an inverse fermion propagator: 

r ^ 2 f i 

G (x, y) 5 a/3 = |i<9 t + — | + \i + v(x) - / d 4 zll (x - z)a(z)^ 5 4 (x - y) 5 a(3 . (5) 

Note that this propagator still depends on the field a but is manifestly diagonal in the 
spin (or flavor) indices. We exponentiate the determinant using detA = exp [Tr In A] and 
reintroduce a source term, now coupled to the a field, to obtain 



Z[J] 



-iE[J] 



T>a exp I gTr In 



G~ (x,y) 



x exp ^- / d 4 x I d 4 y cr(x)U (x — y)a{y) s j exp (i / d 4 xJ{x)a{x) j 



(6) 



Below we will use a generating functional with an external source J coupled to an auxiliary field rather 
than to ip^ip. These functionals generate the same observables; see Ref. |3j| for a discussion of the 
relationship between them. 
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where the spin/flavor trace has been performed in the first term, leading to the spin/flavor 
degeneracy factor g (e.g., for electrons or neutrons, g = 2 and for symmetric nuclear matter, 
g = 4). The remaining trace in the first line of Eq. is over space-time. 

We define the "classical field" cr c (x) in the presence of J(x) by the ground state expecta- 
tion value of a(x): 



a c (x) = (Q\a(x)\Q)j 



SJ(x) 



d E[J} = -t T ^\nZ[J] 



5J(x) 



and define the effective action 

T[a c ] = -E[J] - [ d 4 xJ{x)a c {x) 



in the usual way ||34j| . The effective action has the property 

-T[a c ) =-J(x), 



5a Jx 



(?) 



(8) 



(9) 



and the solutions for J(x) = represent the stable quantum states [34]. At the minimum 
a® (with J(x) = 0) of a uniform system, the energy density S of the ground state is related 
to the effective action by 

T K]\ J=0 = -VTS, (10) 

where VT is the space-time volume of the system. More generally, at finite density we must 
examine spatially dependent cr° to find the absolute ground state. 

To evaluate T[a c ] in a loop expansion, we write a = cr c + r] and expand Eq. @ in quantum 
fluctuations rj around the classical field a c . Once again we seek a Gaussian path integral, 
this time in terms of r), while using a source coupled to rj to treat the residual //-dependent 
terms perturbatively (by removing them from the path integral in favor of functional deriva- 
tives with respect to the source). The expansion, integration, and subsequent Legendre 
transformation are standard and we simply quote the result for the effective action from 
Refs. g|,|3§: 

T[a c ] = ^■Ti\n[Gj I 1 (x,y)} + ^ J d 4 x j ' d 4 y a c (x)U (x - y)a c (y) 

+ -Trln [D~ 1 (x,y)j + / d 4 x5C[o- c ] + (connected IPI-diagrams) . (11) 



In Eq. flTT|) , we have introduced the inverse Hartree propagator (which depends on a c but 
not rf) 



^ 2 



d 4 zUo(x — z)a c (z) 5 4 {x — y) , 



(12) 



G H 1 {x,y) : 

and the inverse a propagator 

D~ 1 (x, y) = -iU (x -y)+g J d 4 z 1 J d 4 z 2 G H {zi, z 2 )U (z 2 - x)G H (z 2 , ^1)^0(^1 - y) , (13) 

which originates with the part of the exponential in Z[J] that is quadratic in rj after ex- 
panding. The counterterm Lagrangian S£ has been included in Eq. ([LI]) for completeness; 
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however, with the regularization/renormalization procedure applied here (DR/MS), we will 
not need it explicitly [[34 



The two propagators are obtained by solving the equations 

j ' d 4 zGjj 1 (x,z)G H (z,y) = 5 4 (x-y), 

Jd 4 zD- 1 (x,z)D a (z,y) = 6\x-y), (14) 

with appropriate boundary conditions for Gu and D a (discussed below). The "connected 
lPI-diagrams" in Eq. (Jllf) are built from the propagator D a for the a field and the vertices 



V m (x 1 ,..,x m ) = d 4 yt... d 4 y m G H (y m ,yi)U (y 1 - x l )..G H {y )U (y ), (15) 



where m = 3, . . . , oo. We follow Ref. [3|[] and use the freedom of a counterterm 5 J to 
ensure that a c is the same to all orders in the expansion (i.e., (rj) = 0), which simplifies the 
bookkeeping. As a consequence, all diagrams that are 1-particle reducible with respect to 
the a propagator are cancelled and only 1-particle irreducible diagrams have to be included 
in calculations of the effective action. In the next section, we show how this expansion, when 
applied to the lowest-order EFT potential, is an expansion in inverse powers of g. 



III. DILUTE FERMI SYSTEMS WITH SHORT-RANGE INTERACTIONS 

In this section, we apply the effective action formalism to the Fermi gas with short- 
range, spin-independent interactions. We calculate the effective action to one loop explicitly 
and demonstrate that the loop expansion of the effective action corresponds to a large N 
expansion. 



A. Short-range EFT and Low-density Expansion 

We consider a general local Lagrangian for a nonrelativistic fermion field that is invariant 
under Galilean, parity, and time-reversal transformations: 



2 



£ = iQ t + J_ ^ °-(^iP) 2 + — {^)\i)W 2 i)) + h.c. 



2M 



2 VT T/ 16 



+ ^(^VV) t -(W) + ... , (16) 

where V is the Galilean invariant derivative and h.c. denotes the Hermitian 

conjugate. The terms proportional to C 2 and C' 2 contribute to s-wave and p-wave scattering, 
respectively, while the dots represent terms with more derivatives and/or more fields. The 
Lagrangian Eq. (fjj]|) represents a particular canonical form, which can be reached via field 
redefinitions. For example, higher-order terms with time derivatives are omitted, as they 
can be eliminated in favor of terms with spatial derivatives. 

To reproduce the results in Ref. 0, we can write a generating functional with the la- 
grangian Eq. (|1^) and Grassmann sources coupled to ^ and ip, respectively. Perturbative 
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P/2 + k P/2 + k' 

X 

P/2-k P/2-k' 





+ V + V + 




k 2 + k' 2 

-i{k'\V EF T\k) -iC -iC^—^r -iC' 2 k-k' 

FIG. 1: Feynman rules for (k'|VEFT|k). P is the total momentum in the center of mass. The spin 
indices have been suppressed. 

expansions for Green's functions follow by taking successive functional derivatives, and the 
ground state energy density follows by applying the linked cluster theorem (see Ref. p5| ). 



The scattering amplitude for fermions in the vacuum is simply related to a sum of Feyn- 
man graphs computed according to this lagrangian. The terms in Eq. (|T6|) involving only four 
fermion fields reduce in momentum space to simple polynomial vertices that are equivalent 
to a momentum expansion of an EFT potential for particle-particle scattering, 

(k'\V EFT \k) = C + C 2 (k' 2 + k 2 )/2 + C^k' • k + . . . , (17) 

(see Fig. [1] for the Feynman rules). Because of Galilean invariance, the interaction depends 
only on the relative momenta k and k' of the incoming and outgoing particles. The co- 
efficients Co, C*2, and C 2 can be obtained from matching the EFT to a more fundamental 
theory or to (at least) three independent pieces of experimental data. It is important to note 
that Eq. ( |TT| ) is not simply the term-by-term momentum-space expansion of an underlying 
potential because the coefficients also contain short- distance contributions from loop graphs. 
By matching to the effective range expansion, we can express the Cn in terms of the effective 
range parameters: 

C = ^ C 2 = C ^ and C> = ^i, (18) 
U M 2 2 M v; 

where a s and r s are the s-wave scattering length and effective range, respectively, and a p is 
the p-wave scattering length. (The implied renormalization prescription here is DR/MS, as 
described in Ref. ||.) 

The energy density at finite density can be calculated from diagrams with no external 
lines using the vertices in Fig. [I| and the propagator 

Co(fco,k)„ a = 5 aS ( «£ - *f >, + - k \ ) , (19) 

\ k - c^ k + it k Q — uj^ — it J 

where = k 2 /(2M). In this approach, which is applicable to a dilute Fermi system at 
T = 0, the finite density boundary conditions are imposed by hand in Eq. (0), rather than 
using a chemical potential and inverting it at the end flOf . The leading contribution to the 



energy density is the kinetic energy of the noninteracting Fermi gas 



d 3 k 3k 2 
E kin = g I j^u^9{k ¥ -k)=n- 1 ^. (20) 
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w CO =- o-o + (T) 

(t) («) 



(») 44 tt => 4 t 4 t + 




(0 («) 

FIG. 2: Hugenholtz diagrams contributing to the energy density of the dilute Fermi gas at O(kp) 
(a) and O(kp) (b). The right-hand side shows the two possible spin-contractions for each of the 
two diagrams. 

where 

n = gk 3 F /(67r 2 ) (21) 

is the number density of the system. In the low-density limit, the leading corrections come 
from the diagrams shown on the left-hand side of Fig. |2|. 

For tadpole diagrams such as in Fig. 0(a), a convergence factor exp(ikor)) must be in- 
cluded, with the limit rj —>■ + being taken after the contour integrals have been carried 
out. This procedure automatically takes into account that such lines must be hole lines. 
Diagram |2|(b) does not need the convergence factor, but has a power law UV divergence 
and must be renormalized. A convenient regulator is dimensional regularization with min- 
imal subtraction (DR/MS) B. For a detailed description of the Feynman rules and the 
calculation of the energy density for a dilute Fermi gas to 0(k F In /cf) ? including three-body 
contributions, see Ref. |J. To O(kp) the energy density is JTU| 2 



S = n 



h ^ 3 + ( g -l){l-(k F a s ) + J_ (H_21n2)(fc F a s ) 2 }] (22) 



2M 



.5 1 37r s 35n 



where the first term is the kinetic energy of the noninteracting Fermi gas and the second 
and third terms are corrections from the diagrams in Fig. 0(a) and (b), respectively. 

In the low-density expansion of Eq. (|22"D , it is assumed that g ~ 0(1), and therefore 
the direct and exchange contributions to Fig. |2|(a) are counted at the same order in the 
kp expansion (and similarly with higher-order contributions). Thus, Hugenholtz diagrams, 
which combine these contributions, are particularly efficient for this expansion. If we study 
systems where g can be considered large, then we will need a new power counting that 
assigns direct and exchange contributions to different orders in the EFT expansion. The 
isolation of different g dependencies is accomplished in the diagrams on the right-hand side 



2 Note that the low-density expansion corresponds to an expansion in powers of and that the overall 
factor of density implicitly contains three powers of fcp. 
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+ 



--o- + 



FIG. 3: Integral equation satisfied by the a propagator (double dashed line). The single dashed 
line indicates the bare a (a single potential exchange) while the solid lines are Hartree propagators. 

of Fig. |2|, which follow from the introduction of an auxiliary field. We consider a consistent 
power counting and summation of such diagrams in the next section. 



B. Effective Action for Short-Range Interactions 



Here we will calculate the energy density for the Fermi gas with short-range interactions 
using the effective action formalism. The loop expansion of the effective action does not 
correspond to a simple low-density expansion but resums an infinite class of diagrams. We 
will only explicitly consider the momentum independent Co interaction from Eq. (|T6|). The 
other C21 vertices and many-body forces can be included perturbatively (see Appendix [A]). 

The interaction Uq from the previous section can then be specified as 



U (x -y) = C S 4 (x - y) 



(23) 

which corresponds to the second term in Eq. (|16"1). 3 Inserting Eq. (|23|) into Eq. (|TlD , we 
obtain 



I> c 
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Trln[G^(x,y)] + 



Cn 



d 4 x a c (xY 



+ -Trln[D; 1 (x,y)] 



+ (connected lPI-diagrams) 
The inverse Hartree and o propagators simplify to 



G H l (x,y) 
D-\x,y) 



id t + 



+ n + v(x) - C cr c (x) 



5 4 {x - y) 



2M 

iC 5 4 (x -y)+ gCl G H (y, x)G H (x, y) 



(24) 

(25) 
(26) 



The er propagator D a (x, y) is given by the sum of all fermion ring diagrams and satisfies the 
integral equation 



(27) 



D a (x, y) = £T^ 4 -y)~ igCo I d 4 z D a (x, z)G H (z, y)G H (y, z) 



3 Strictly speaking there is a contribution proportional to ip'ip from normal ordering, but this can be 
absorbed into a counterterm for the chemical potential. 



10 



(a) (b) 

FIG. 4: Vertices (a) V3 and (b) V4 for the lPI-diagrams from Eq. (p8|). The solid lines indicate 
Hartree propagators Gh while the double dashed lines indicate the a propagator D a . 



which is illustrated in Fig. [| The solid lines denote Hartree propagators Gn{x,y) while the 
double dashed line denotes the a propagator D a (x, y). The single-dashed lines represent the 
bare a propagator given by the inhomogeneous term in Eq. (p7|); it corresponds a single a 
exchange. The bare a propagator carries a factor 1/Cq while each coupling of the bare a 
to fermions carries a factor Co, such that a four-fermion interaction is proportional to Co as 
required. 

The "connected lPI-diagrams" are built from D a (x,y) and the vertices 



V m (x l ,..,x m ) = (C ) m G H ( 

•£"mi • • • Gj[ yX m .. \ . X m ) 



(28) 



where m — 3, . . . , 00. In Fig. f|, we show the lowest-m vertices, V3 and V4, as an example. 
The solid lines indicate Hartree propagators Gh while the double dashed line indicate the 
a propagator D a . The "connected lPI-diagrams" can be arranged in a consistent loop 
expansion, where the corresponding expansion parameter is the total number of loops minus 
the number of fermion loops [see Eq. (|50|)1. The leading corrections to the terms explicitly 
given in Eq.([24j) start at two loops. As an example, we show the "connected lPI-diagrams" 
at two loops in Fig. |5|. All other diagrams have at least three loops. 

The status of the loop expansion as an expansion in 1/g can be seen by returning to the 
generating functional of Eq. (fj), but now specialized to the contact interaction of Eq. (f23|): 



Z[J] = e ~ iE[J] = J Vo exp \ #Trln 



^ 2 

i9t + 2M +fl ~ C ° a ^ 



x exp ^ i I d 4 x ( -Coa(x) 2 + J(x)a(x) 



(29) 



If we scale Co and a according to Co = co/g and o = ga', then we find that all of the g 
dependence reduces to a single overall factor multiplying the path integral exponent. Thus, 
a saddle point evaluation of the integral is just the loop expansion with the loop counting 
parameter equal to 1/g (which plays the same role as h in the usual discussion) p5| . Below 
we will examine the contributions order by order. 



C. Leading-order Effective Potential 

In the following we carry out the effective action formalism for a uniform Fermi gas with 
short-range interactions. For simplicity, we work without an explicit chemical potential and 
incorporate the appropriate finite-density boundary conditions in the propagator ||, [L0| . We 
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(a) 




+ 




(&) 



FIG. 5: "Connected lPI-diagrams" contributing to the effective action at two loops. The solid 
lines indicate Hartree propagators Gh while the double dashed lines indicate the a propagator D a . 



also set the external potential v(x) to zero, since we are interested at present in an infinite 
system. We assume that the ground state is uniform. As a consequence, the classical field 
a c (x) is independent of x and we can set a c (x) = a c in Eq. (f25|). We reassess this assumption 
in Sect. [TV]. 

The first step is to calculate the contribution to the effective potential given by the first 
line in Eq. (|24] ) . We refer to this contribution as the tree- level contribution since it does not 
depend on D a . The TrlnD" 1 in the second line of Eq. ([24]) gives the one-loop part of the 
effective potential. Since Gjj 1 is diagonal in momentum space, the Trln can be evaluated as 



TrlnG^ 1 



(30) 



where 

e p = p 2 /(2M) + C a c (31) 

and VT is the spacetime volume. The dpo integral can be regularized by taking a derivative 
with respect to e p and integrating back after the dpo integral is carried out. (The lower 
integration limit can be taken to be independent of p and will not contribute to the d 3 p 
integration in dimensional regularization.) The resulting integral has the form of a tadpole 
diagram with a modified propagator. The appropriate boundary conditions can be included 
by the replacement 



1 



Po 



G H (p) 



6(p-k F ) 6(k F -p) 



Po~e p + it po 



e p - it 



(32) 



Applying the convergence factor exp(i^p ) an d taking the limit 77 — > + after the contour 
integral has been carried out, we obtain 



-TrlnG^ 1 

i H 



-gVT 



d 3 p 



p_ 

2M 



C a c \8(k F -p). 



(33) 



The d 3 p integral is now immediate and the full tree contribution to the effective potential 
from the first line of Eq. (123) is 



3 kp 



5 2M 



n — Cqg c n 



(34) 
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By requiring the effective action to be stationary, 



5a, 



(J,. 



n . 



(35) 



and substituting cr° back into Eq. 



we obtain the energy density at tree level 



£ 



to 



'3 k 2 
n I - 



5 2M 



n- 



2M 



3 2 , 

5 3tc 



(36) 



Note that this expression includes the Hartree term [diagram (i) in Fig. 0(a)] while the 
Fock term [diagram (ii) in Fig. 0(a)] will appear in the Trln of the inverse a propagator 
in the second line of Eq. (|2"3|). [Diagram (i) of Fig. |2|(b) will also appear in the Trln at 
one-loop order while diagram (ii) of Fig. ||(b) will appear at two loops, as seen from the first 
equality of Eq. (|50f).1 As a consequence, the loop expansion of the effective potential does 
not correspond to the low-density expansion of the dilute Fermi gas where the Hartree and 
Fock terms appear at the same order |§. 

Instead, the loop expansion of the effective potential corresponds to a 1/N expansion 
where N is the spin degeneracy factor g. This can be seen as follows: We take the limit 
g — > 00 for the energy per particle E/N = £/n by defining a new coupling 



Co = g C , 

and keep Co and kp fixed as g — > 00. Using Eqs. (|36"1 , |37D , we find 



E 

N 



3 kp 



+ 



Co 



LO 



5 2M 2 6tt 



2 • 



(37) 



(38) 



which is 0(1) in the 1/g expansion. Note that Eq. (37) implies that the scattering length a s 
is 0(1/ g) in the large g limit. Thus, scattering in free space is perturbative, since a particle 
can scatter off only one other particle, while the many-body problem is nonperturbative as 
g — » 00 since a particle can scatter off g other particles. 



D. Next-to- leading Order Effective Potential 

We proceed to calculate the one-loop contribution to the effective potential given by 
Trhi-D" 1 in the second line of Eq. ([24]) . We will show below that this term is of 0(1/ g) 
and constitutes the first correction to Eq. (|36"D . After performing a four- dimensional Fourier 
transform 

Gh(x x ,x 2 ) = J ^G H (p)e-^-^ , (39) 
the contribution from the second line in Eq. (|24|) to the effective potential becomes 

FnloWc] = 5lW>[/ 0ji e ' i9<Xl ' X2) (-<Co+0C« j -^G H (p + q)G H (p) ) j 

= \VT J (0i ln (-^o + ^C o 2 J J^G B (p + q)G H (p)) , (40) 



13 



where Gh{p) is the Hartree propagator with the appropriate boundary conditions given in 
Eq. ©. 



It is customary to define the polarization insertion [ID 



/d f d 

j^iGoip + q)G (p) = -ig J J^Gh(p + q)G H (p) , (41) 

where Go(p) is given in Eq. fllPP and the second equality is valid for a uniform system, for 
which a c is a constant that can be eliminated by a shift in pq. By taking a derivative with 
respect to Co in Eq. ([40|) and integrating back, we obtain 



r \„ l 1 vt f C ° dy f dAq yUo{q) (A9) 

2 Jo V J (2vr) 4 l-?/n (g) 

where we have used that Jd 4 q/(2n) A = in dimensional regularization. The first two 
terms in the expansion of 1/(1 — yUo(q)) in powers of yUo(q) correspond to the Fock term 
(diagram (ii) in Fig. |](a)) and diagram (i) in Fig. y(b). While all higher-order terms in the 
expansion are UV finite and can be summed straightforwardly, those two terms are special: 
The diagrams in Fig. |2|(a) are tadpoles that require a convergence factor and the diagrams 
in Fig. H|(b) have power law UV divergences and need to be renormalized. Therefore, we 
subtract those two terms from Eq. (0) and calculate them separately using the methods of 
Ref. [1. 

The remaining finite terms are easily summed by calculating 

f \« 1 1 VT [ C ° dy f d ' q ^ n ^)) 3 (a* 
Tnlo[Vc\ = --VT / — / — — — — . (43) 



2 Jo y J (27r) 4 l- y n (g) 
Performing the dy integral in Eq. (|43"D leads to 



TnloW] = ~VT J 



f A ' ' ^ ' ' ' - (44) 



(27T) 



CoTlo(q) + ^(C n (g)) 2 + ln(l - C n (g)) 



and it is straightforward to verify that the first two terms in the integrand simply cancel 
the first two terms in an expansion of the logarithm. Note the similarity of Eq. fl44T) to the 



expression [Eq. (12.53)] for the correlation energy of a degenerate electron gas in Ref. [10 



Because Tnlo does not depend on a c (for uniform background fields), it does not have to 
be minimized and is, up to a factor of — VT, equal to the contribution to the energy density 
at 0(1/ g). Since the energy density is real, we can write 



poo poo 

£nlo = ~ 777 / v 2 dv I dv, 



4tt 3 M 
+ arctan 



C Imn (vo,^) + C lmIlo(vo,v)ReIlo(vo,v) 
C lmUo(vo,v) 



(45) 



C ReTlo(vo,v) - 1, 

where v = |v| and (fo,v) are dimensionless variables related to the four-momentum (qo,q) 
via 

= |q|/&F and v = Mqo/kp. (46) 
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Adding the contributions of the Fock term and diagram 
to the energy density is given by 



of Fig. ||(b), the full contribution 



£ 



NLO 



n 



'2M 



-—(k F a s )+g—(ll-2]n2)(k F a s 



NLO 



(47) 



Using Eqs. (|18|) and (|37|), one can verify that the contributions to the energy per particle 
from Eq. ( (4"T|) are indeed of 0(1/ g). To this order, similar RPA contributions appear in the 
modified loop expansion discussed in Refs. and We will calculate Snlo numerically 
and discuss the energy density below. In the next section we comment on the higher-order 
contributions. 



E. Higher-order Effective Potential 

The higher-order contributions come from connected lPI-diagrams built from the vertices 
in Eq. ( ffijp and the a propagator 

= £(i + *q, l^ H{p + q) G M y^± T - 7 l m . (48) 

The diagrams occuring at two-loops are shown in Fig. |5[ Using Eqs. ( |37|) and (fP|) it is easy 
to see that D a ~ O(g). It is then straightforward to verify that the lPI-diagrams shown in 
Fig. |5] are of 0(1/ g 2 ) in the energy per particle. In fact, these are the only contributions at 
this order, with all other diagrams being suppressed by at least one more power of 1/g. A 
general argument can be made as follows: 4 Fermion loops and a propagators contribute a 
positive power of g while a a-fermion-fermion vertex is proportional to Co and contributes 
a negative power. The contribution to the energy per particle of diagram with Lp fermion 
loops, V a-fermion-fermion vertices, and I a a propagators scales as g v with 

v = L F - V + I a - 1 . (49) 

We define the number of a loops as the total number of loops minus the number of fermion 
loops L a = L — Lp. The number of a propagators must be half the number of a-fermion- 
fermion vertices l a = V/2, and twice the total number of internal lines is equal to three 
times the number of a-fermion-fermion vertices, 21 = 3V. Using the general topological 
relation relating the total number of loops L, vertices V, and internal lines I of a given 
diagram: L = I — V + 1, we can rewrite Eq. ( fl9D as 

v = Lp — I a — 1 = —(L — Lp) = —L a , (50) 

which proves that the loop-expansion of the effective action is equivalent to a 1/g expansion 
of the energy per particle if no many-body forces are present. 

An important question is at what order in 1/g three-body forces enter. In Appendix 
[A], we discuss how many-body forces can be included in our formalism. Here we use a 
renormalization argument to show that they will contribute at 0(l/g 2 ). In Fig. [| we 



4 The argument does not apply in the present form if many-body forces are present. However, we will use 
a renormalization argument below to identify at which order three-body forces contribute. 
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FIG. 6: Diagrams with logarithmically divergent contributions to the effective action at two loops. 
The solid lines indicate Hartree propagators Gh, the double dashed lines indicate the full a propa- 
gator D a , and the single dashed lines indicate the bare a propagator (a single potential exchange) . 



illustrate two particular contributions that are contained in the two-loop diagrams from 
Fig. [5] by expanding the o propagator. These two diagrams have logarithmic UV divergences. 
This can be seen most easily by cutting the hole lines to obtain the corresponding diagram 
in free space. The resulting diagrams are 1-particle irreducible and describe three-particle 
scattering. They have the same UV divergences as the finite-density diagrams in Fig. ||. 
Proper renormalization of this UV divergence requires a three-body contact interaction of 
the type (ip^ip) 3 . This was first noticed for a system of bosons in Ref. [|1^] and later applied 
to fermions in Ref. ||. For the 1/g expansion to be consistent, the three-body counterterm 
must appear at the same order as the diagrams it renormalizes, so three-body forces have 
to enter at 0(1/ g 2 ). This is manifest in the renormalization-group equation for the running 
of the non-derivative contact three-body force D (see Ref. 0), which depends on Cq. 



IV. ENERGY DENSITY AND STABILITY ANALYSIS 



A. Calculation of Snlo 



We now turn to the calculation of Snlo- In order to evaluate the integral in Eq. fljSp, we 
need the polarization insertion LTo. The calculation of the polarization insertion is standard 
and can be found, for example, in Ref. HlOfl . Here we only quote the final result: 



ReIIo(fo,f) 



Imn (t>o, v) 



9MkF h + ±-(l-(v /v-v/2f)ln 



4tt 2 



--(l-( Vo / v + v /2) 2 ) In 



1 + (v Q /v - v/2) 



1 - (vq/v - v/2) 
(v /v + v/2) 



(v /v + v/2) 



gMfa 



iirv 



(1 - (v /v - v/2) 2 ) ; 



v 2 V 2 
v > 2, — + v>v > — 



gMk F . 2 v 
-— (1 - {v v - v/2) ) ; 

67TV 

gMk F 



2Vn 



inv 



v < 2, < v < v 



V V 

v < 2, v H > v > v 

v 2 



(51) 
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FIG. 7: The universal function H{x) for —4 < x < 4. The solid line shows the result of the 
numerical calculation. The inset shows H(x) for — 1 < x < 1 in more detail. 



= ; otherwise . (52) 

For the calculation of Snlo it is convenient to define the dimensionless variable 

_ gMkp g 
x = - C = -k F a s , (53) 

47T 2 71 

and the universal function 

F (v ,v) = ^^U (v ,v). (54) 



The integral in Eq. ( ]45[) can then be written as 

h 5 

4tt 3 M 



£nlo = . o F . , / v 2 dv I dv n 



x ImF (t>o, v) + x 2 Im F (v , v)Re F (v , v) 
x Im Fq(v , v) 



— arctan 



7.5 
F 



4tt 3 M 



1 — x Re Fo(uo, f ) 

i?(x) (55) 



where H(x) is a universal function of x alone. Due to the complicated form of Uo{vo,v), 
Eq. ( [5*5] ) has not been evaluated analytically. Since £nlo is manifestly finite, however, H(x) 
can be calculated numerically. In Fig. 0, we show the results of a numerical evaluation of 
H(x). The inset shows H(x) in the perturbative region — 1 < x < 1 in more detail. The 
local minimum for negative x allows for a self-bound uniform system. (Repulsion at NLO 
was also observed in the modified loop expansion of Ref. Pj3.) 
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0.0 
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Ms 



FIG. 8: Energy per particle in units of 1/Ma 2 s at leading order (LO) and next-to-leading-order 
(NLO) in the 1/g expansion, for g = 4 and g = 10. 



The energy per particle through NLO can be written in terms of g and x as 

E _ kl IY3 2 \ 1 (3 „, s 2.4 
N ~ 2M 



- + -x ) + - ( -H(x) - ~x+ — (11 - 21n2)x 2 
5 3 J g \7T 3 35 v 7 



(56) 



In Fig. D, we show the energy per particle in units of the Fermi energy for two different 
values of g: g = 4 and g — 10. The dashed lines show the LO results while the solid lines 
show the NLO order results. The difference between LO and NLO gives an indication of 
where the expansion breaks down. The NLO result oscillates around the LO result in the 
vicinity of kpa s = and grows without limit as kpa s — > ±oo. As g increases, the minimum 
becomes shallower and moves toward smaller values of kpa s . 

An important consideration is whether the energy calculated above actually corresponds 
to the true ground state of the system. There are two issues here: whether the state that 
evolves adiabatically from the noninteracting ground state is unstable to pairing near the 
Fermi surface and whether a homogeneous system is stable. There is no pairing instability if 
we restrict ourselves to repulsive interactions (e.g., a positive, natural-size scattering length), 
but pairing is inevitable with an attractive interaction (e.g., a negative scattering length or 
x < 0). We will investigate the role of pairing in the EFT framework in future work. We 
note here that more conventional investigations typically find that the impact of pairing on 
the bulk properties (on which we focus) is small for natural scattering lengths [[IT| . 

In order to investigate the stability of the homogeneous state, we calculate the compress- 
ibility || 



d ( 2 d E 



A ' = M"W' (57) 

When the compressibility is positive, the system is stable with respect to infinitesimal long- 
range fluctuations in the density. When the compressibility is negative, however, long-range 
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fluctuations will grow exponentially. As a consequence, the system will cluster and the 
homogeneous state is not the ground state. Using Eq. (p6|), the compressibility at leading 
order in the 1/g expansion is 



K 



LO 



2kp ( gkpa s 



3M 



7T 



which turns negative for 



x < x 



crit 



(58) 



(59) 



As a consequence, the homogeneous system is unstable if x < — 1/2. This constraint was 
previously obtained in Ref. [ 5jJ by identifying the radius of convergence for the exact fermion 
RPA energy, as well as from a variational argument. At next-to- leading order in the 1/g ex- 
pansion, the compressibility K^lo follows from Eqs. (p6|) and (p7|) but we have not obtained 
an analytic expression. Furthermore, the point where the compressibility becomes negative 
now depends on g and k F a s separately. Nevertheless, the numerical results for K NLO are 
close to the leading order values fl59l) . For example, the next-to-leading order compressibility 
turns negative if kpa s < —0.73,-0.37,-0.15 for g = 2,4,10, respectively, as compared to 
facts < —0.79, —0.39, —0.16 at leading order. This means that the self-bound uniform state 
given by the minimum of E/N in Fig. ||| is not a stable ground state. 

The above constraints can also be obtained directly from the effective action by looking for 
minima in Eq. (|24|) corresponding to a nonuniform state (with a c = a c (x)). The saddle points 
of the effective action determine where such minima appear. This analysis also determines 
the density fluctuation modes that destabilize the homogeneous state. In leading order, we 
require: 







5a c (x)5a c (y) 



FloWc(z)] 



5a c (x)Sa c (y) 



9 -^HG H \x 1 y)} + < ^ 



d xaJx) 



where a c is a function of x. Using the identity 

5 2 



: Trln[G^] 



-C%G H {x,y)G H {y,x) 



, (60) 



(61) 



5a c (x)5a c (y) 

performing a Fourier transform with respect to x — y, and setting q = 0, we obtain the 
equation 



= l-C II (g = 0,g) = l + :r<h 



kh 



4Af 



In 



q/{2h 



1 + 9/(2**) 



(62) 



whose solutions determine the long-range density fluctuation modes that destabilize the 
homogeneous system. The constraint obtained from solving Eq. ( |B"2"D agrees with Eq. (pP[). At 
this order, the stability analysis in the effective action formalism is equivalent to identifying 
the poles in the o propagator D a (q = 0, q) from Eq. (f48l), as is done in conventional many- 
body approaches [10 , 25] . In principle, it is straightforward to extend the stability analysis 
to next-to-leading order by requiring 







d 2 



Sa c (x)5a c (y) 



{^loWc(z)} + r NLO [cr c (z)}} . 



(63) 



However, this analysis involves the evaluation of two-loop diagrams and is beyond the scope 
of this work. 
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V. BOSE LIMIT 



The EFT expansion is not limited, in general, to small values of x. For example, Eq. (p6[) 
implies that the 1/g expansion is valid for large x if g is also sufficiently large. The Bose 
limit, in which x —>■ oo, provides an instructive example. 

The energy density of a Bose system can be obtained from the energy density of a Fermi 
system by taking the limit g — ► oo, k F — ► 0, but keeping n = gk F /(6Tr 2 ) constant P5[ . Here 
g is the degeneracy of an artificial "flavor" quantum number. This can be understood intu- 
itively starting from a noninteracting Fermi system with large degeneracy. If the degeneracy 
g is greater than the number of particles, the spatial state will be a symmetric wave function 
with all particles in the lowest momentum state. This is the same as the wave function of 
the Bose system with the same number of particles. If the interacting state evolves adia- 
batically from the noninteracting state, then the limiting Fermi system is the same as the 
interacting Bose system times a totally antisymmetric "flavor" wave function that has no 
physical consequences [[32 . 



The Bose limit can be taken directly in Eq. (0), leading to the mean-field energy for a 
dilute Bose system. All terms vanish except for the part of the second term proportional to 
g, and the result is 

d 2ira s n 2 , , 

— ■ (64) 

in agreement with the literature [TIJ. For higher-order terms, the analysis is more com- 
plicated since some contributions diverge as g — > oo and therefore need to be resummed. 
However, this is precisely the same resummation that is implemented in the one-loop effective 
action from the previous section. 

In the Bose limit, the polarization insertion Hq [see Eq. (£!)] has a simple analytic 



expression. In order to evaluate Eq. fl£2]) most easily, we take the Bose limit of Ho before 
performing the dy and d A q integrals. The polarization insertion can be written as 

n (go,q) = g / -—,6(k F - P ) — ■ — ■ . (65) 



(2tt) 3 \q + lu p - lu p+ci + ie q - u p + u; p _ q - ie 

If we let |q| = q ^> k F (or equivalently take k F 0), the 9 functions involving q can be 
dropped. Expanding 

^ P - ^p±q = -<*> q (l + 0(p/ q)) , (66) 
the polarization insertion in the Bose limit becomes [Q 

n B (g ,q) = 9 [ i%0(h-p) ( ^— -) 

J (2vr) d V<?o-^q + 2e qo+Uq-ieJ 

2U « n (67) 



(g - ^ q + «e)(<?o + ^ q - ie) ' 

Thus Hq has simple poles in q given by the noninteracting single-particle kinetic energy. 

The ring sum leading to the integrand in Eq. ( f42"[) shifts the poles to the Bogoliubov 
quasiparticle energies: 

yUo(q) V 2uqny 2u; q m/ 



1 - yn (g) n \q) -y ql-ul- 2u (l ny + ie ql~e\ + ie J 
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where 

e q = \Ju\ + 2w q ra/ . (69) 
The dqo integral in Eq. (j42|) is simply evaluated by contour integration and leads to 



n b _ J_nlo _ n [ c ° f d 3 q o; q _ n f Co f d\ q 2 
1 " VT ~ 2 J dy J (2vr)3 eq " 2 J * V J (2nf q ^ q 2 + AM ny ' ( j 

Note that an expansion of the final integrand in powers of y would generate infrared (IR) 
divergences. The resummation of the one-loop fermion effective action is therefore required 
for bosons even in the low-density limit. This summation of ring diagrams because of IR 
divergences is analogous to the calculation of the correlation energy for a uniform electron 



gas, only in that case it is the interaction that is the source of the IR divergences |L0 . 

The d 3 q integral in Eq. (|70D can be evaluated in dimensional regularization with minimal 
subtraction. In particular, the formula 

d d k e r(-d/2)r(i/2 + d/2) /2 



{2n) d kVk 2 + A 2 (47r)<V 2 r(d/2)r(l/2) 
can be derived and applied for d = 3 |34| . The end result is 



n , r ,,/o^5/2 2na fl n 2 128 



^^M^^^^, (72) 



in agreement with the well-known result [ T(| . 



In principle, our numerical evaluation of the universal function H(x), illustrated in Fig. [7], 
should match on to the Bose limit for large, positive x > 10. We have not been able to 
demonstrate this matching. Indeed, the numerical continuation to large x of H(x) exhibits 
an entirely different behavior (e.g., it reaches a maximum and then turns negative, going 
to minus infinity with a different asymptotic power of x then predicted by the Bose limit). 
It is not clear whether the problem is one of numerics (our calculations x ^ 10 show severe 
round-off errors) or a more fundamental problem. Further investigations to resolve this issue 
are in progress. 



VI. SUMMARY 

In this paper, a system of fermions with short-range interactions at finite density is studied 
using the framework of effective field theory. The effective action formalism for fermions with 
auxiliary fields leads to a loop expansion in which particle-hole bubbles are resummed to all 
orders. For spin-independent interactions, the loop expansion is equivalent to a systematic 
expansion in 1/N, where "iV" is the spin-isospin degeneracy g. Thus we have a double 



5 If one starts from the subtracted expression (^3|) instead of Eq. (f42|), the same result is obtained with a 
direct integration. The separately calculated terms proportional to fcpa s and (k-pa s ) 2 in Eq. ( |4^ ) vanish 
in the Bose limit, while the subtractions implicit in Eq. ( f43| ) cancel linear and cubic UV divergences. If 
dimensional regularization with minimal subtraction is used as above, these UV divergences are subtracted 
automatically and one can use Eq. ( fl2|) directly. 
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expansion, in 1/g as well as kpa s . The loop expansion differs from the dilute expansion 
described in Ref . , even at leading order. At next-to-leading order, the expansion requires 
Hartree plus RPA contributions to the ground state 0, [35], [4(J. 

The formalism enables us to examine uniform systems for potential self-bound solutions. 
At next-to-leading order in 1/g, the nonperturbative resummation leads for the uniform 
system to a non-trivial and non- analytic dependence on the Fermi momentum k-p. There is 
a self-bound minimum within the radius of convergence of the EFT for sufficiently large g, 
but a stability analysis reveals that it is unstable to density fluctuations (e.g., clustering). 

An interesting limit of the large g expansion takes g to infinity and kp to zero, with the 
density n held fixed. This limit reproduces the energy expansion for a dilute Bose gas, with 
the non-analytic dependence on \fn generated by the infinite summation of ring diagrams. 
We have not succeeded in reproducing the Bose limit starting from our explicit calculations 
for fermions, which may simply reflect the numerical difficulties in carrying out the limit 
but might also mean that there are subtleties we have not recognized. Work is in progress 
to resolve this issue. 

While the present investigation illustrates some features of a systematic nonperturbative 
treatment of a Fermi system, we need to extend our treatment to include pairing and large 
scattering lengths to make contact with the physical systems of greatest interest. Investiga- 
tions in these areas are in progress as well as work to adapt the effective action formalism 
to finite systems in the form of density functional theory [[O . 
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APPENDIX A: HIGHER-ORDER TWO-BODY AND MANY-BODY FORCES 

In this appendix, we generalize the discussion of the 1/N expansion to consider higher- 
order two-body interactions (C2, C f 2 , and so on) as well as many-body forces (D and so 
on) that are present in a general EFT Lagrangian for the dilute system (see Ref. || for 
definitions and details): 

(^) t (V'V 2 V)+ h.c. 

+ ■ (VWO - ^(V>W + • ■ • • (Al) 

8 

There are various different ways of incorporating these interactions into the loop expansion, 
which correspond to different reorganizations of the original perturbative expansion, with 
different resummations of diagrams. The details of the system under consideration determine 
the power counting of the higher-order contributions; the most appropriate resummation will 
implement that power counting. 

At the simplest level, we can incorporate all higher-order vertices perturbatively. This 
means that the only infinite summations based on large g will be of the same subdiagrams 
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FIG. 9: (a) New diagrams and (b) a propagator and vertex corrections generated by the pertur- 
bative inclusion of higher-order terms in the effective Lagrangian (|A1|) . 



with Cq vertices we have already considered. We return to the generating functional in 
Eq. @, but now with the full EFT Lagrangian for a dilute system with short-range inter- 
actions, and with Grassmann sources £t and £ coupled to the fermion fields: 

Z[^,£\ = e ~ iE ^'® = [ V^Vi!e^ d4x[c+ ^ (x) ^ x)+ ^ {x) ^ x)] , (A2) 



where £ is given by Eq. fllEp. [Spin indices are implicit in Eq. ( JA2[ ).] Functional derivatives 
with respect to the Grassmann sources can be used to generate the n-point fermion Green's 
functions. If we isolate the quadratic part of the Lagrangian, this procedure generates the 
perturbative expansion of these functions in terms of the noninteracting fermion propagator 
(at finite chemical potential). The sum of closed, connected diagrams generated this way 
reproduces the perturbative expansion of the energy density from Ref. §|. 

In the present context, we use the Grassmann sources to remove all interaction terms from 
the path integral Ci(^, ip) — > £/(— 2<5/<5£, —iS/Stf) except for the Co contact interaction. To 
carry out the procedure in Sect. II leading from Eq. (H) to Eq. @, we need to incorporate 
the new Grassmann source terms. This is straightforward, since the identity 

^ G -^ + ^ + £ty = (^t + tfOG- 1 ^ + GO - , (A3) 

and a shift in ip and ^ leaves the same integral as before (leading to the same Tr In term) 
plus a new term: 



exp 



i / d A xd A y£(x)G(x,y)£(y) 



(A4) 



The new term still depends on a through G(x,y). 

As in Eq. ([□]), we expand a about a c in Eq. (|A4T), which has the effect of expanding 
G about Gh- The Grassmann sources £ and ^ in Eq. (|A4j) have to match up with the 
derivatives in Ci(—i8/8£, —i6/6$) because the sources are set to zero in the end. This dic- 
tates how the higher-order vertices (e.g., Ci and D ) appear in Feynman diagrams attached 



to Gh propagators. The quadratic term in rj can be directly incorporated into Eq. ([26 ) 
while the higher-order terms generate additional connected lPI-diagrams. (For consistency, 
one should expand the a propagators perturbatively in the higher-order interactions.) The 
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77- independent part of Eq. (|A4|) , which has G — > Gh, generates tadpole-like diagrams with 
the new vertices, such as those illustrated in Fig. ^|(a), as well as higher-order contributions. 
The new diagrams based on the rj expansion can be constructed by breaking every possible 
Gh line (including those buried in the a propagator) and connecting the broken ends to the 
new vertices. Examples of cx-propagator and V3 vertex corrections are shown in Fig. [|(b). 

An alternative approach to incorporating higher-order two-body and many-body inter- 
actions is to generalize the a Lagrangian in Eq. fl2D| ) to include a 3 and higher-order self- 
interactions, as well as terms with gradients (but not time derivatives) acting on a. If 
we put back the fermion path integral and eliminate the a field by iteratively applying its 
Euler-Lagrange equation via field redefinitions, we reproduce an infinite subset of the terms 
in the general dilute effective Lagrangian. For example, we recover the Do vertex and the 
combination of C2 and C 2 corresponding to a term proportional to V^fy) • V(ip'ip). In 
general, we find all interaction terms with ip^ipa to some power and gradients acting on such 
factors (these are terms that depend only on the momentum transfer between the fermions). 

The generalized a Lagrangian can be expanded about a c as before, with the new terms 
treated to all orders (analogous to the treatment of the linear sigma model in Ref. ||34|| ) or as 
perturbative corrections. In the former case, the strict power counting in which all diagrams 
contribute at a given order in 1/g is violated unless the counting of the coefficient of the 
a 3 term, D oc D , is promoted to 0{l/g 2 ) from 0(l/g 4 ) in the perturbative case. The 
latter case is similar to the discussion above, except that a given higher-order term in the 
generalized a Lagrangian will be equivalent to a selective summation of higher-order terms 
expressed in the basis of fermion fields alone. 
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